Measurements of the observed cross sections for e+e exclusive light hadrons 
containing Kg meson at s/s = 3.773 and 3.650 GeV 
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By analyzing the data sets of 17.3 pb~^ taken at ^/s — 3.773 GeV and of 6.5 pb~^ taken 
at = 3.650 GeV with the BES-II detector at the BEPC collider, we measure the observed 
cross sections for the exclusive light hadron final states of KgK~ tv^ , KgK~ n^-K° , KgK~Tv^-K^-K~ , 
KgK~ Tv^Tv^TT~Tv'^ , KgK~Tv^-K^-K'^Tv~TT~ and KgK~TT'^Tv^Tv'^ produced in e"'"e~ annihilation at the 
two energy points. We set the upper limits on the observed cross sections and the branching fractions 
for ^(3770) decay to these final states at 90% C.L.. 



I. INTRODUCTION 



The -0(3770) resonance is expected to decay almost 
entirely into DD meson pairs since its width is almost 
two orders of magnitude larger than that of ■0(3686) 
In recent years, the study of the 0(3770) non-DD de- 



nium energy region due to the existing puzzle that about 
38% of ■0(3770) does not decay into DD meson pairs 
0- To understand the possible excess of the 0;(377O) 
cross section relative to the DD cross section, BES and 
CLEO Collaborations made many efforts to study the 
0^(3770) non-DD decays. The CLEO Collaboration mea- 
sured the e+e" ^0(3770) non-DD cross section to be 
(-0.01±0.08+S ai) nb 11. While the BES Collaboration 



2 



measured the branching fraction for tjj(3770) non—DD 
decay to be (15 ± 5)% [1, S, H, S which indicates 
that, contrary to what is generally expected, the ip{3770) 
might substantially decay into non—DD final states or 
there are some new structure or physics effects which may 
partially be responsible for the largely measured non-DD 
branching fraction of the -0(3770) decays [1, BES 
Collaboration observed the first non—DD decay mode for 
■!/'(3770) —^ J/ipTr^'K~ , and measured its decay branching 
fraction to be S[V'(3770) ^ J/V-tt+tt-] = (0.34 ± 0.14 ± 
0.09)%Jll[l|. This was confirmed by CLEO Collabo- 
ration [19|. Latter, CLEO Collaboration observed more 
0(3770) exclusive non-DD decays, 0(3770) ^ J/ipn°TT°, 
J/^n°, J/^ri [il,7Xcj(J = 0,l,2) [H, HI and 0r; 0, 
etc. Summing over these measured branching fractions 
yields the sum of the branching fractions for the 0(3770) 
exclusive non-DD decays not more than 2%. In ad- 
dition, BES and CLEO Collaborations also attempted 
to search for other 0(3770) exclusive charmless decays 
[13, m, m, [13, mi tll, [H, [H. However, the existing 
results can not clarify the possible excess. For better un- 
derstanding the origin of the possible excess, search for 
more -0(3770) exclusive charmless decays will be helpful. 

In this Letter, we report measurements of the ob- 
served cross sections for the exclusive light hadron final 
states of KgK~TT'^ (Throughout the Letter, charge con- 
jugation is implied), i^^i^T^Tr+Tr", KgK~TT~^TT^TT^ , 

X^X^TT+Tr+TT^Tr*', KgK^TT^TT^TT^TT^TT"' and 

KgK'TT^Tr'^Tr'^ at the center-of-mass energies of 3.773 
and 3.650 GeV with the same method as the one used 
in our previous works [l^ [13, [IH . With the measured 
cross sections at the two energy points, we set the upper 
limits on the observed cross sections and the branching 
fractions for 0(3770) decay to these final states. The 
measurements are made by analyzing the data set of 
17.3 pb-i collected at = 3.773 GeV [called as the 
0(3770) resonance data] and the data set of 6.5 pb~^ 
collected at y/s = 3.650 GeV (called as the continuum 
data) with the BESII detector at the BEPC coUider. 



II. BESII DETECTOR 



The BES-II is a conventional cylindrical rnagnetic de 
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tector that is described in detail in Refs 
12-layer vertex chamber (VC) surrounding the beryllium 
beam pipe provides input to the event trigger, as well as 
coordinate information. A forty-layer main drift chamber 
(MDC) located just outside the VC yields precise mea- 
surements of charged particle trajectories with a solid an- 
gle coverage of 85% of 47r; it also provides ionization en- 
ergy loss (dE/dx) measurements which are used for parti- 



cle identification. Momentum resolution of 1.7%a/1 -I- 
{p in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering electrons are obtained for the data taken at 
-yi = 3.773 GeV. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOP) 
of charered narticles with a resolution of about 180 ns for 



electrons. Outside the TOP, a 12 radiation length, lead- 
gas barrel shower counter (BSC), operating in limited 
streamer mode, measures the energies of electrons and 
photons over 80% of the total solid angle with an energy 
resolution of (Te/E = 0.22/a/B {E in GeV) and spatial 
resolutions of — 7.9 mrad and — 2.3 cm for elec- 
trons. A solenoidal magnet outside the BSC provides a 
0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return and serve to identify muons with 
momentum greater than 500 MeV/c. They cover 68% of 
the total solid angle. 



III. 



EVENT SELECTION 



In the 



reconstruction 



KlK- 

thc Kl and 



KlK- 



and K^K- 



of the 



K%K--K 



final states, 



decays of 



mesons are reconstructed through the 
7r+7r~ and tt" — > 77. 
Por each candidate event, we require that at least four 
charged tracks are well reconstructed in the MDC with 
good helix fits, and the polar angle of each charged track 
satisfies \cos9\ < 0.85. The charged tracks (except for 
the Kg meson reconstruction) are required to originate 
from the interaction region Vxy < 2.0 cm {Vxy < 8.0 
cm) and \Vz\ < 20.0 cm, where Vxy and 1141 are the 
closest approaches in the xy-plane and the z direction, 
respectively. 

The charged particles are identified by using the dE/dx 
and TOP measurements, with which the combined con- 
fidence levels CLjr and CLx for pion and kaon hypothe- 
ses are calculated. The pion and kaon candidates are 
required to satisfy CLjr > 0.001 and CLk > CLy^, re- 
spectively. To reconstruct Kg mesons, we require that 
the TT^TT^ meson pairs must originate from a secondary 
vertex which is displaced from the event vertex at least 
by 4 mm in the xy-plane. 

The photons are selected with the BSC measurements. 
The good photon candidates are required to satisfy the 
following criteria: the energy deposited in the BSC is 
greater than 50 MeV, the electromagnetic shower starts 
in the first 5 readout layers, the angle between the photon 
and the nearest charged track is greater than 22° [2^.[27|. 
and the opening angle between the cluster development 
direction and the photon emission direction is less than 

37° [li,[13. 

Por each candidate event, there may be several differ- 
ent charged and/or neutral track combinations satisfying 
the above selection criteria for exclusive light hadron fi- 
nal states. Each combination is subjected to an energy- 
momentum conservation kinematic fit. Por the processes 
containing 7r° meson in the final states, an additional con- 
straint kinematic fit is imposed on tt'^ — 77- Candidates 
with a fit probability larger than 1% are accepted. If 
more than one combination satisfies the selection criteria 
in an event, onlv the combination with the loneest decav 
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distance of the reconstructed Kg mesons is retained. 

To suppress the background from the DD decays, we 
use the double tag method [1^ to remove the DD events. 
For example, for the KgK~TT~^TT'^ final state, we exclude 
the all possible events from DD decays by rejecting those 
in which the D and D mesons can be reconstructed in 
the decay modes of D~ KgK^ and D~^ 7r"'"7r°, 
D- ^ i4r-7rO_and D+ K'^n+ , D° K'^tt° and 
— > K^TT^ [28]. For the other final states, the events 
from DD decays are suppressed similarly. The remaining 
contaminations from DD decays due to particle misiden- 
tification or missing photon(s) are accounted by using 
Monte Carlo simulation, as discussed in Section fVl 



IV. DATA ANALYSIS 
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In the data analysis, these processes containing Kg 

meson in the final state are studied by examining 

the invariant mass spectra of the tt+tt" combinations 

satisfying the above selection criteria for the Kg meson 

reconstruction. The invariant masses of the tt+tt" 

combinations are calculated with the momentum vectors 

from the Kg reconstruction. Figure [1] shows the resulting 

combi- 
+ 



distribution of the invariant masses of the tt^tt 
nations from the selected candidates for the K^K 



KlK- 



In 





s 

K^K~ TT^ ■K'^ -K^ TT~ TT~ and K^K'TT^-K^TT^ final states 
each figure, the peak around the K^ nominal mass indi- 
cates the production of e^e~ — > exclusive light hadrons 
containing K'g meson. Fitting the tt+tt^ invariant mass 
spectra with a Gaussian function for the A'^ signal and 
a fiat background yields the number of the events for 
each process observed from the V'lSTTO) resonance data 
and the continuum data. In the fit, the K^ mass and 
its mass resolution are fixed at the values obtained by 
analyzing Monte Carlo samples. 



V. BACKGROUND SUBTRACTION 

Some other events may contribute to the selected can- 
didate events for e+e" ^ / (/ represents exclusive light 
hadron final state). These include the events from J/ij] 
and V'(3686) decays due to ISR returns, the events from 
the other final states due to misidentifying a pion as a 
kaon or reverse, and the events from DD decays. The 
number of these contaminations should be subtracted 
from the number N°^^ of the candidates for e+e^ f. 
The estimation of them can be done based on Monte 
Carlo simulation. The details about the background sub- 
traction have been described in Ref. [l^ . For each back- 
ground channel except DD decays, 50,000 or 100,000 
Monte Carlo events are used in the background estima- 
tion. The Monte Carlo sample of each different back- 
ground channel is from ten to several thousands times 
larffer than the data in size. 



Fig. 1: The ■k'^tt invariant mass spectra of the candidates 
for the (a) K°sK-Tr+, (b) K^K--jt+tt°, (c) K'^R- tt+ty+tv- , 
(d) K'^K-Tr+TT+-!v-TT°, (e) K'^K-n+Tr+Tr+n-Tr- and (f) 
A'sA'^tt+tt'V final states selected from the i/)(3770) reso- 
nance data (left) and the continuum data (right). 



Monte Carlo study shows that the contaminations from 
V'(3770) ^ J/tPn+TT-, ^^(3770) ^ J/tpn°7r°, ^{3770) -> 
J/ijn^ and ■0(3770) jxa {J = 0,1,2) can be ne- 
glected. 

Even though we have removed the main contamina- 
tions from DD decays in the previous event selection (see 
section HlH) ■ there are still some events from DD decays 
satisfying the selection criteria for the light hadron final 
states due to particle misidentification or missing pho- 
ton(s). The number of these contaminations from DD 
decays are further removed by analyzing a Monte Carlo 
sample which is about forty times larger than the V'(3770) 
resonance data. The Monte Carlo events are generated 
as e+e~ DD at ^/s = 3.773 GeV, where the D and D 
mesons are set to decay into all possible final states with 
the branching fractions quoted from PDG 

Subtracting the number of these contaminations 
from the number N°^^ of the candidate events, we ob- 
tain the net number A^"''' of the signal events for each 
process. For the K^K'tt'^ and Ar^A'~7r+7r+7r+7r~7r~ fi- 
nal states, for which only a few signal events are observed 
from the continuum data, we set the upper limits A^^p on 
the number of the signal events at 90% C.L.. Here, we 
use the Feldman- Cousins method ^] and assume that 
the background is absent. The numbers of N°^^, 
and iV"°' (or N^p) are summarized in the second, third 
and fourth columns of Tabs. [TlandlTTl For each process, 
the background events in the ■)/'(3770) resonance data are 
dominant by DD decays and -0(3686) decays. While, 
there is no DD decay in the continuum data, and the 
ip{3Q8Q) production cross section at y/s — 3.650 GeV is 
much less than that at = 3.773 GeV. So. the num- 
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ber of the background events in the continnum data can 
almost be neghgible. 



VI. 



RESULTS 



Monte Carlo efficiency 



To estimate the detection efficiency e for e 



f, we 

use a phase space generator inchiding initial state radi- 
ation and vacuum polarization corrections [30] with 1/s 
ener gy d ependence in cross section. Final state radia- 
tion [3l| decreases the detection efficiency not more than 
0.5%. Detailed analysis based on Monte Carlo simulation 
for the BES-II detector [s^l gives the detection efficiencies 
for each process at -^i = 3.773 and 3.650 GeV, which are 
summarized in the fifth columns of Tabs. Ill and [Til where 
the detection efficiencies do not include the branching 
fractions for Kg tt+tt" and 7r° 77, B{Kg tt+tt") 
and B{tt^ 77). 



states at y/s = 3.650 GeV are set with Eq. ([T]) by sub- 
stituting A^"''* with iV"P/(l - Asys), where TV"? is the 
upper limit on the number of the signal event, and Agys 
is the systematic error in the cross section measurement. 
Inserting the corresponding numbers in the equation, we 
obtain the upper limits on the observed cross sections for 
e+e^ — !■ KgK^TT'^ and e~^e~ —> KgK~TT~^'!T'^'!T~^TT~Tr~ 
at — 3.650 GeV, which are also listed in Tab. HH 



C. Upper limits on the observed cross sections and 
the branching fractions for 7/)(3770) f 

If we ignore the possible interference effects between 
the continuum and resonance amplitudes, and the dif- 
ference of the vacuum polarization corrections at y/s = 
3.773 and 3.650 GeV, we can determine the observed 
cross section cr^(377o)— / for ■0(3770) — > / at \/s = 3.773 
GeV by comparing the observed cross sections a^:^'^'^'^}^ 



3.650 GoV 



for e+e f measured at -^/s = 3.773 



and a 

and 3.650 GeV, respectively. It can be given by 



B. Observed cross sections 



cr, 



0(3770)-/ 



3.773 GcV 



fo 



„ ^3.650 GoV 

X CTe+e-^f i 



(2) 



Let Bt^o = S(7r" ^ 77) for the modes of 
K'^K-TT+Tr° and K^K-tt+tt+tt-tt^ , B^o = B^{tt^ 
77) for the mode of Kc;K~tt^tt^tt'^ and B^^o = 1 for 
K'^K-TT+, K'^K-^TT+n+TT- and K^K-tt+tt+tt+tt-tt- , 
where B{Tr'^ — s- 77) is the branching fraction for the decay 
of tt" — > 77, then the observed cross section for e~^e~ f 
can be determined by 



N" 



X BtjO ' 



(1) 



where £ is the integrated luminosity of the data set, A^"°* 
is the number of the signal events, e is the detection ef- 
ficiency and B{Kg — > tt+tt") is the branching fraction 
for the decay of Kg — s- tt+tt". Inserting these numbers 
in Eq. ([T]), we obtain the observed cross sections for 
each process at y/s = 3.773 and 3.650 GeV. They are 
summarized in Tabs. [T] and |TT1 where the first error is 
statistical and the second systematic. In the measure- 
ments of the observed cross sections, the systematic er- 
rors arise from the uncertainties in integrated luminos- 
ity of the data set (2.1% 0, Q), photon selection (2.0% 
per photon), tracking efficiency (2.0% per track), parti- 
cle identification (0.5% per pion or kaon), kinematic fit 
(1.5%), Kg reconstruction (1.1% [27|). branching frac- 
tions quoted from PDG ^ (0.03% for B{tt° 77) 
and 0.07% for B{K'^ tt+tt")), Monte Carlo modehng 
(6.0% [llSIIH), Monte Carlo statistics (1.4% ~ 4.4%), 
background subtraction (0.0% ~ 3.0%) and fit to mass 
spectrum (0.4% ~ 8.5%). Adding these uncertainties in 
quadrature yields the total systematic error Agys for each 
mode at ^/s — 3.773 and 3.650 GeV. 

The upper limit cr"+^_ , , on the observed cross sec- 
tions for the KiK-Tr+ 



and K°cK- 



where fco — 3.650^/3.773^ is the normalization factor to 
consider the 1/s cross section dependence. The results 
are summarized in the second column of Tab. IIIIl where 
the first error is the statistical, the second is the inde- 
pendent systematic arising from the uncertainties in the 
Monte Carlo statistics, in the fit to the mass spectrum 
and in the background subtraction, and the third is the 
common systematic error arising from the other uncer- 
tainties as discussed in the subsection B. 

The upper limit on the observed cross section 
'^V'(3770)-./ ^-(3770) ^ / at = 3.773 GeV is set 
by shifting the cross section by 1.64(7, where a is the to- 
tal error of the measured cross section. The results on 



(3770)-/ 



are summarized in the third column of Tab. 



The upper limit on the branching fraction S"'^ 



V(3770)-/ 

/ is set by dividing its upper limit on 
y by the observed 

(7.15 ± 0.27 ± 0.27) nb [H for 

the -0(3770) production at y/s = 3.773 GeV and a factor 
(1 — Aa°^f^^^jQ~^) , where Acr^'^g^-^g^ is the relative error of 



for V(3770) 

the observed cross section t7j^^377o) 
cross section o-'^(^y70) 



the i-r°^^ 

" V(3770) ■ 



The results on B^^^^^^^^^j are summarized 



in the last column of Tab. IIIIl 



VII. SUMMARY 

In this Letter, we present the measurements 
of the observed cross sections 



X^AT-tt+ttO, 



A:gA:"7r+7r+7r- 



for K°gK-T:+, 
K^K-n+n+TT-n^. 



final 



Ar^A'^7r+7r+7r+7r~7r" and ATgAT^Tr+Tr^Tr" produced in 
e+e^ annihilation at y/s — 3.773 and 3.650 GeV. These 
cross sections are obtained bv analvziner the data sets of 
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17.3 pb"^ taken at = 3.773 GeV and of 6.5 pb'^ at 
Vi = 3.650 GeV with the BES-II detector at the BEPC 
coUider. By comparing the observed cross sections for 
each process measured at = 3.773 and 3.650 GeV, 
we set the upper hmits on the observed cross sections 
and the branching fractions for V(3770) decay to these 
final states at 90% C.L.. These measurements provide 
helpful information to understand the mechanism of the 
continuum light hadron production and the discrepancy 
between the observed cross sections for DD and ?/'(3770) 
production. 
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Tab. I: The observed cross sections for e^e" — > exclusive light hadrons at ^/s = 3.773 GeV, where A*'"'''' is the number of events 
observed from the 1/1(3770) resonance data, A'^'' is the total number of background events, A''"'^* is the number of the signal 
events, e is the detection efficiency, Aaya is the relative systematic error of the observed cross section, and a is the observed 
cross section. 





^obs 






e{%) 




<7°'==[pb] 


K%K- 






18.4 ±4.6 


0.1 ±0.0 


18.3 ±4.6 


10.02 ±0.14 


10.7 


15.2 ± 3.8 ± 1.6 


K%K- 




^° 




41.2 ±6.6 


1.1 ±0.2 


40.1 ±6.6 


3.52 ±0.08 


11.6 


96.2 ± 15.9 ± 11.1 


K%K- 






TV 


40.0 ±6.5 


1.0 ±0.2 


38.9 ±6.5 


3.56 ±0.06 


14.2 


91.5 ± 15.3 ± 13.0 


K%K- 




7r+ 




24.5 ±5.2 


1.5 ±0.3 


23.0 ±5.2 


0.77 ±0.03 


15.2 


253.0 ±57.1 ±38.4 


K%K- 








4.8 ±2.2 


0.3 ±0.1 


4.5 ±2.2 


0.84 ±0.03 


18.4 


44.4 ±21.9 ±8.2 


K%K- 


7T^ 






19.8 ±4.9 


2.8 ±0.5 


17.0 ±4.9 


0.99 ±0.04 


14.3 


147.0 ±42.4 ±21.0 



Tab. II: The observed cross sections for e^e^ exclusive light hadrons at ^/s = 3.650 GeV, where A^"*"* is the number of 
events observed from the continuum data, A'"^ is the upper limit on the number of the signal events, a"'' is the upper limit on 
the observed cross section set at 90% C.L., and the definitions of the other symbols are the same as those in Tab. H] 





yyobs 


ATb 


Af"'='(or Af"P) 


e{%) 




cr°'^=(cr"P)[pb] 


K°sK- 


^^ 




2 


0.0 


< 5.91 


10.55 ±0.15 


12.9 


< 14.3 


K%K~ 


7^^ 


TT 


7.7 ±2.9 


0.0 


7.7 ±2.9 


3.62 ±0.09 


11.6 


47.9 ±18.0 ±5.6 


K%K^ 


7^^ 


- + - 
TT TT 


13.4 ±3.8 


0.0 


13.4 ± 3.8 


3.66 ±0.06 


14.1 


81.4 ±23.1 ±11.5 


K%K- 


7^^ 


- + - 
TT TT TT 


4.6 ±2.5 


0.0 


4.6 ±2.5 


0.87 ±0.03 


17.2 


119.0 ±64.7 ±20.5 


K°sK- 


^^ 


- + + — - 
TT TT TT TT 





0.0 


< 2.44 


0.95 ±0.03 


18.1 


< 69.7 


K%K- 


7^^ 


TT TT 


3.3 ±2.0 


0.0 


3.3 ±2.0 


1.12 ±0.05 


14.1 


67.1 ±40.7 ±9.5 



Tab. Ill: The upper limits on the observed cross section ct^p^tto)-,/ \fs — 3.773 GeV and the branching fraction ^3^^3779)^^ 
for ?/)(3770) — > / are set at 90% C.L.. The 0-^(3770)-./ is calculated with Eq. ([2]), where the first error is the statistical, the 
second is the independent systematic, and the third is the common systematic error. The upper * denotes that we neglect the 
contributions from the continuum production. 



Decay Mode 


O"V{3770)^/ (Pb) 


'^^(3770)-./(Pb) 


cup 

'-'i/.(3770)^/ 




15.2 ±3.8 ±0.2 ±1.6* 


< 22.0 


3.2 X 10"^ 


K°sK-7T+n'' 


51.4 ±23.2 ±2.6 ±5.8 


< 90.7 


13.3 X 10"^ 


K%K-Ti+iT+Ti- 


15. 3 ±26. 5 ±2. 2 ±2.1 


< 59.0 


8.7 X 10"^ 


K%K-Ti+n+Ti-n° 


141.6 ±83.2 ±14.7 ±20.7 


< 284.3 


41.8 X 10"^ 


K%K- TT+ 71 + 11+ TT-TT- 


44.4 ±21.9 ±2.1 ±7.9* 


< 82.7 


12.2 X 10"^ 


K%K-^+iv'\'' 


84.2 ± 57.2 ± 8.3 ± 11.2 


< 180.4 


26.5 X 10"^ 



